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Abstract The structure of electric current and magnetic helicity in the solar 
corona is closely linked to solar activity over the 11-year cycle, yet is poorly 
understood. As an alternative to traditional current-free "potential field" ex- 
trapolations, we investigate a model for the global coronal magnetic field which 
is non-potential and time-dependent, following the build-up and transport of 
magnetic helicity due to flux emergence and large-scale photospheric motions. 
This helicity concentrates into twisted magnetic flux ropes, which may lose equi- 
librium and be ejected. Here, we consider how the magnetic structure predicted 
by this model — in particular the flux ropes — varies over the solar activity cycle, 
based on photospheric input data from six periods of cycle 23. The number of 
flux ropes doubles from minimum to maximum, following the total length of 
photospheric polarity inversion lines. However, the number of flux rope ejections 
increases by a factor of eight, following the emergence rate of active regions. 
This is broadly consistent with the observed cycle modulation of coronal mass 
ejections, although the actual rate of ejections in the simulation is about a 
fifth of the rate of observed events. The model predicts that, even at minimum, 
differential rotation will produce sheared, non-potential, magnetic structure at 
all latitudes. 

Keywords: Magnetic fields. Corona; Magnetic fields. Models; Solar cycle. Mod- 
els; Coronal Mass Ejections, Theory 



1. Introduction 

The number of sunspots has long been known to vary with a period of about 
11 years. Since the early 20th century it has been recognised that sunspots 
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carry magnetic flux originating from the Sun's interior. Once this magnetic 
flux emerges into the corona, it dominates both the structure and evolution 
of the plasma there. This is evidenced by 11 -year variation s in many phenom- 
ena, including solar flares ( Aschwanden , 20051). prominences ( d'Azambuia , 1955 



Hansen and Hansenl. 19751). X-ray flux ( Pevtsov and Actonl. 2001 ), corona l mass 



ejections fCM Es: IWebb and Howardl Il994t [Cr emadcs and St. Cvil l2007l ). and 
the solar wind (jRichardson. Wang, and Paularenai. .20011) . 

Despite the connection to observed phenomena, the coronal magnetic field 
cannot usually be measured directly d ue to the tenuous nature of the plasma, so 
only i solated measurements exist {e.g 



Lin. Penn, 



inuous nature oi tne plasm 
and Tomczvd . 120001: ICasir 



et g^H 



20031 : iTomczvk et all |2007|) . Usually, it must be extrapolated from the routinely 



observed magnetic field in the solar photosphere, perhaps using images of high- 
temperature loops to constrain the field topology. This t opology is often found to 

conta in non-zero electric currents {e.g. in X-ray sigmoids; Canfield. Hudson, and McKenzi^ .[ 
19991) . and extrapolation of suc h magnetic fields is not yet robust, even when 



limited to a single active region ( DeRosa et al. . 20091 ). On the global scale, which 
concerns us in this paper, it is usual to assume a current-free magnetic field (a 
potential field). This is uniquely determined in the corona given an observed 
radial component on the photosphere and vanishing hor izontal components at an 



uppe r "source surface" (usually placed at r = 2.5Rr^: lAltschuler and Newkirk 



19691) . Even the more realistic models that solve for MHD equilibria in the 



corona usually assume a potential magnetic field which is then per turbed to be 



consistent with a given t hermodynamic structure and solar wind (jRilev et al 



20061 : ICohen et all 120071 ). Whilst the potential field model has been successful 
in describing certain large-scal e aspects of coronal magnetic struc ture — such as 
the locations of coronal holes ( Wang. Hawlev. and Sheelev . 1996[) — it does not 
include the development of currents and free magnetic energy that are needed 
to model solar eruptions. In fact, the highly s heared magnetic fields observed in 
long- lived filament channels all over the Sun ( Martin. Bilimoria. and Tracadad . 



19941 ) show that, even outside active regions, the assumption of a potential field 



can be inadequate. 

As a first step toward understanding the structure of currents and magnetic 
helicity in the global corona, we have investigated an alternative model for 
the coronal magnetic field, based also on observational input of the photo- 
spheric magnetic field. In this model, the large-scale mean magnetic field in 
the corona evolves in time through a quasi-static relax ation, in re sponse t o flux 
emergence and to surface motions Ivan Ballegooiien, Priest, and M acka^ 



Yeates. Mackav. and van Ballegooiieni 120081 : lYeates and Mackavi l2009bf ). The 



consequence is that, unlike in a sequence of potential field extrapolations, current 
and magnetic helicity are generated in and transported through the corona. 
Flux cancellation above photospheric polarity inversion lines (PILs) then leads 
to the concentration of helicity in either sheare d arcades or twisted magnetic 
flux ropes (through the mechanism described bv lvan Ballegooiien and Martens . 
19891) . Comparison of the simulated magnetic field direction at these locations 



with observations of filament chirality has demonstrated that, despite its sim- 
plicity, the model is able to reproduce the general structure of the magn etic field 
at most of these locations ( Yeates. Mackav. and van Ballegooiieni 2008). 
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In this paper, we use the quasi-static model to simulate the coronal mag- 
netic field during six distinct periods over cycle 23. Our aim is to consider 
how the magnetic field structure predicted by this model varies over the 11- 
year solar activity cycle. We focus on a particular aspect of the field structure: 
the formation and ejection of magnetic fiux ropes. These are not only observed 
in the real corona, bu t their ejection is suggested to give rise to CMEs {e.g., 
iGibson and Fan . 2006 ). By their very nature, flux ropes cannot be modelled 
by potential field extrapolations, so their locations and ejections have not been 
considered in previous models of the global corona. The formation of fiux ropes 
in the present quasi-static model has been describe d for a simple two-bipole 
configuration bv lMackav and van Balleeooiien J 2006h. and for t he global corona 
over a particular 5-month period by Yeates and Mackav ( 2009b[ ). Here we extend 
the latter study to simulate different phases of the solar cycle. 



2. Coronal Magnetic Field Simulations 

Our numerical simulati ons of the global coronal magnetic fi e ld ar e based on 
the mean-field model of van Balleeooiien. Priest, and Mackavl (|2000f ). and were 



described in detail by lYeates and Mackav ( 2009bf ). The kev features are as 
follows: 

1. The large-scale mean-field in the corona evolves through the non-ideal MHD 
induction equation, with a turbulent diffusivity to parametrize the effect of 
small-scale fluctuations. 

2. We impose an artificial "magneto- frictional relaxation" velocity so that the 
coronal field responds to photospheric driving by evolving through a sequence 
of force-free equilibria. 

3. The photospheric dr iver is th e stan dard surface flux transport model for the 



radial magnetic field ISheelevI ()2005() 



4. New active regions are inserted in the form of idealised magnetic bipoles with 
properties (location, size, tilt, and magnetic flux) based on observations. The 
3D bipoles are given a non-zero helicity. 

The simulations use a stretched spherical coordinate grid between radii Rp, 

and2.5i?0 with an effective resolution of 1° at the equator (see Yeates. Mackav. and van Ballegooiienl [ 



20081 ). As observational input we use normal-component synoptic magnetograms 
from the US National Solar Observatory at Kitt Peak. These are used in two 
ways: to generate the initial condition for each run (a potential fie ld extrapola- 



tion) , and to determine the properties of newly-emerged bipoles (see Yeates. Mackav. and van Balleeooiien^ 



(j2007[ ) for details). In this paper, we consider six simulation runs, each of 5.5 
Carrington rotations in duration. These periods are labelled A to F, and shown 
in Table [TJ They are chosen based on the availability of the Kitt Peak data, so 
as to cover a range of phases of solar cycle 23 from one solar minimum to the 
next. The fourth column of Table [T] shows the total number of bipoles emerged 
in each simulation run, while the right-most column indicates the instrument 
used for the magnetogram observations; the older vacuum telescope (KPVT) 
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Table 1. Simulation periods. | 



Period 


Carrington rotations 


Dates 


Bipolos 


Instrument 


A (minimum) 


1911.5 to 1917 


12-Jul-96 to 04-Jan-97 


16 


KPVT 


B (rising phase) 


1948.5 to 1954 


17-Apr-99 to ll-Oct-99 


117 


KPVT 


C (maximum) 


1962.5 to 1968 


03-May-OO to 27-Oct-OO 


122 


KPVT 


D (maximum) 


1974.5 to 1980 


26-Mar-Ol to 19-Sep-Ol 


109 


KPVT 


E (declining phase) 


2018.5 to 2024 


08-Jul-04 to 02-Jan-05 


52 


SOLIS 


F (minimum) 


2067.5 to 2073 


06-Mar-08 to 30-Aug-08 


14 


SOLIS 



was replaced in 2003 by SOLIS (Synoptic Optical Long-term Investigations of 
the Sun). 

Figure [1] shows example synoptic magnetograms from each of the six periods 
A to F (left column) , along with snapshots of each simulation run (middle col- 
umn). These snapshots are taken on the days corresponding to 180° Carrington 
longitude in the synoptic magnetograms. 

An important feature of our non-potential model is the ability to include non- 
zero magnetic hclicity in emerging bipoles, through a parameter /3 in the math- 



ematical expression used to model them (jYeates. Mackav. and van Ballegooiien 



(j2008[ ). Equations (6)-(9)). The magnitude and, especially, the sign of this emerg- 
ing bipol e helicity have been shown to influence both the chirality of filament 
channels (Yeates. Mackav. and van Ballegooiien, 2008; Yeates and Mackav. 2009al) [ 



and the ejection rate of magnetic flux ropes ( Yeates and Mackav . 2009bl) . Unfor 



tunately, the optimum value of the helicity parameter for each bipole is poorly 
constrained by the available observations, and must be arbitrarily chosen. In 
this paper, we depart from our previous practice of using the same value of /3 for 
all bipoles in each hemisphere. Instead, we give each bipole a randomly chosen 
P value. This is motivated by the limited existing observations of active region 
helicity, which show variation in magnitude and sign — both within individual 
active regions and between different regions ( Hagino and Sakurail, 2004) — and 



also a negative gradient of helicity with latitude (jPevtsov. Canfleld. and Metcall 



19951 ) ■ To approximate these observations, we select the /3 value for a bipole at 
latitude A° from a normal distribution with mean /3o = — 0.4A°/25° and standard 
deviation 0.4. The selected values of /3 for each bipole arc shown as a function 
of latitude in the right column of Figure [TJ for each simulation period. Here the 
dashed line shows the mean value /3o(A°) of the normal distribution. Note that, 
in these plots, the sizes of the symbols indicate the relative amounts of magnetic 
flux in each bipole. Visible also is the sunspot "butterfly diagram" whereby 
emerging active regions progress toward the equator as the cycle progresses. 
While the random choice of helicity in each emerging bipole clearly limits the 
application of the model to the detailed structure of particular active regions, 
we are interested here in the global magnetic picture, and how this varies over 
the solar cycle in this model. 



SOLA: arycycle.tex; 25 March 2010; 0:57; p. 4 



Solar Cycle Variation of Magnetic Flux Ropes 



CRi917 








* 










CR1954 




it M 


















MMlMMiMMiMi 




m -0.5 







0.0 
m -0.5 : 
-1.0 



"~ V ^ 




» ^ 

Period A ^ 








Period B 









** * 

Period C 








Period D )K 






* 


Period E x 




«^ 




Period F 







Latitude (deg) 



Figure 1. Summary of the six coronal magnetic field simulations. Left column shows example 
synoptic magnctograms from NSO/Kitt Peak in each period, middle column shows snapshots 
of the corresponding simulations, and right column shows hclicity of emerging bipoles as a 
function of latitude. In the middle column, shading shows polarity of radial solar surface 
magnetic field, and lines show selected coronal field lines. In the right column, sizes of asterisks 
are proportional to the magnetic flux of each region. 
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Figure 2. Cycle variation of simulated magnetic field properties: (a) properties of the pho- 
tospheric magnetic field, and (b) properties of the 3D coronal magnetic field. Mean values of 
each quantity are taken over the last fOO days of each simulation period, and error bars give 
standard deviation between different days. Each quantity is normalized by its mean value for 
period A. Grey shading shows monthly smoothed sunspot number from SIDC. 



3. Solar Cycle Variation 

3.1. Global Magnetic Field Properties 

3.1.1. Magnetic Flux 

It is clear from the magnetograms in Figure [T] that there is more magnetic flux 
through the photosphere at solar maximum (periods C and D) than at solar 
minimum (periods A and F). In our simulations, the total (unsigned) radial 
magnetic flux through the photospheric boundary increases by a factor of 4 
between minimum and maximum, from 1.4 x 10^'^ Mx in period A to 5.5 x 10^"^ Mx 
in period C. This relative increase is shown by the solid line in Figure [Ifa). 
Here we plot the mean value over the final 100 days of each simulation period, 
normalised by the value for period A. 

3.1.2. Complexity 

Not only is there more photospheric flux at maximum, but it is divided between 
more positive and negative regions, owing to the larger number of active regions 
emerging. In the simulations, the latter corresponds to the number of bipoles 
inserted, shown in the fourth column of Table [T] and by the dotted line in Figure 
mja). Again, the values are normalized by that for period A (16 bipoles). The 
rate of bipole emergence increases by a factor between 7 and 8 from minimum 
to maximum, similar to the observed sunspot number (shown by grey shading). 

An alternative measure of the magnetic complexity is the total length of 
photospheric PILs at any one time in the simulation. This quantity increases 
only by a factor of about 1.5 from minimum to maximum, as shown by the 
dashed line in Figure [2ja). Of course, the PIL length depends on the spatial res- 
olution of the magnetic field. Here we consider PILs of the large-scale mean field 
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van Ballegooiien. Priest, and Mackav , 2000() , which remain well-defined entities 



over many days and are observed to have physical relevance as the locations 
above which filaments form. 

3.1.3. Non-Potentiality 

The deviation of our simulated magnetic field from potential may be charac- 
terised by integrating either magnetic energy or current over the 3D simulation 
volume. The total parallel current |jo • B0/-B0I (at a single time) increases by 
a factor of between 5 and 6 from minimum to maximum. The mean values for 
100 days of each simulation period are shown by the red curve in red in Figure 
[21Jb). The additional coronal currents at solar maximum originate ultimately 
from the greater bipole emergence rate. Emerging bipolcs inject current both 
directly, because they emerge twisted in our simulation, and indirectly, through 
their interaction with nearb y regions and the further sh earing of their fields by 
subsequent surface motions (jYeates and Mackav . l2009bl) . 



The total magnetic energy Eup is shown in blue in Figure [H^b). As expected, 
it increases from minimum to maximum, from 8.9 x 10"^^ ergs in period A to 
6.7 X 10'^'^ ergs in period D — a factor of about 8, a similar relative increase to the 
bipole emergence rate. If we perform a sequence of potential-field source surface 
(PFSS) extrapolations from the simulated photospheric field, we may calculate 
the corresponding PFSS energy Ep at each time (shown in green in Figure ^jp) . 
The free energy of the simulated non-potential field may then be computed by 
subtracting Ep from -Enp- Shown by the triangles and dot-dashed line in Figure 
[2Kb), this free energy Eup — Ep increases by a factor 15 from minimum to 
maximum, a greater relative increase than any other quantity. The actual mean 
values of free energy are 1.3 x lO'^^ ergs in period A and 2.0 x lO'^'^ ergs in period 
D. 

Since the PFSS energy also increases over the cycle, a relative measure of the 
non-potentiality of the field is obtained by normalising the free energy by the 
PFSS energy. The resulting "relative free energy" {Eup — Ep)/Ep is shown by 
the dashed line in Figure ^b). Here the actual value for period A is 0.18. Unlike 
the original free energy, the relative free energy does not show a clear modulation 
with solar activity. So although the total magnetic energy in the corona is much 
lower at minimum, the field in our simulation remains globally non-potential. 
To illustrate this point. Figure [3] compares the magnetic field on the last day 
of period A, for the non-potential simulation (left) and a PFSS extrapolation 
from the same radial photospheric field (right). Here the Sun is viewed from the 
Southern hemisphere. Coloured lines show coronal magnetic field lines crossing 
various PILs, traced from the same locations on the photosphere in each field. It 
is clear that the field across most PILs is much more sheared in the non-potential 
simulation than in the PFSS extrapolation, even though the field strength is 
weak during this minimum period. This demonstrates how magnetic helicity is 
still generated in the corona at minimum, and not merely by emergence in the 
few remaining active regions. In the simulations, shearing of remnant fields by 
differential rotation generates helicity on a large scale, as evidenced by the south 
polar crown in Figure |31 
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Non-potential Simulation PFSS Extrapolation 




Figure 3. Contrast between non-potential simulation and potential field source-surface ex- 
trapolation at solar minimum (period A). The photosphcric radial magnetic field is the same 
in both cases, shown by contours (white shading and solid lines for positive, gray shading and 
dashed linos for negative). Coloured lines show selected coronal magnetic field lines traced 
from the same photosphcric starting points in each case; closed field lines are in blue (lighter 
shades for field lines reaching higher) and open field lines arc in orange. 



3.2. Magnetic Flux Ropes 

Our method for identifying magnetic flux rope s in simulated magne t ic field data 
is a slight modification of that described by lYeates and Mackav ( 2009b ). As 



explained in that paper, we initially identify points on the computational grid 
satisfying certain criteria — so-called flux rope points — before using a clustering 
algorithm to group these points into flux ropes. 

To identify flux rope points, we compute the vertical magnetic tension force 
Tr — (Bo • V)i3or/Mo and pressure gradient Pr ~ — 9rSo/(2/io) at a subset of 
points in the computational domain between r — Rq and r = 1.44i?0, and select 
points by the following five criteria: 

P,(r,_i,0„(^,;) < -O.4Bo2(r,_i,0„0,), (1) 

Pr{r,+ u0,,(f,,) > OAB^{n+u0^,(|>^), (2) 

Tr{n^i,e„(f>,) > O.4B2(r,_i,0„0,), (3) 
Tr{n+u0,,(j),) < -OABl{n+i,e,,(f>,), (4) 
|jo-Bo| > a*Bl (5) 

where a* = 0.7 x fO~^m^^. The factor Bq has been introduced in criteria ([T]) 
to (H]) to make them independent of magnetic field strength, since the original 
criteria were developed for a period of high solar activity and found to miss flux 
ropes with very weak field at solar minimum. Th e new criteria are calibrate d 
so as to reproduce the earlier flux rope statistics of lYeates and Mackav (l2009bl) . 



That paper covers the same period as period B in this paper, but the simulation 
in this paper differs because the emerging bipoles have been given random twists. 
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Figure 4. Cycle variation of simulated magnetic flux ropes and observed CMEs: (a) normal- 
ized by values for period A, and (b) actual ejection rates per day. Values from simulations 
are averaged over the last 100 days of each period (error bars give the standard deviation 
between diff'erent days). Observed CME rates from the CDAW and CACTus catalogues are 
shown averaged over 100-day periods. Estimation of their errors is described in the text. Grey 
shading shows the monthly smoothed sunspot number. 



rather than aU having the same value of j3. This results in slightly fewer flux rope 
points overall and slightly fewer eruptions, due to the lower net helicity. 

Comparing the six simulations in this paper, the number of flux rope points is 
found to vary over the solar cycle. The mean number of flux rope points per day 
for each simulation period is shown in FigurelUJa) by the dashed line, normalized 
by the mean value for period A (692 points). It increases by a factor of about 2.5 
from minimum to maximum. The dotted line in Figure EJl^a) shows the number 
of flux ropes present at any on e time, determined on each day by clustering the 
flux rope points (as described Yeates and Mackav . 2009bl ). Normalised by the 
mean value for period A (19.2 ropes), this shows a similar factor-two increase. 

This doubling of the number of flux ropes present at any one time may be 
understood from the increase in total photospheric PIL length (which was by a 
factor of about 1.5). This is illustrated by Figure[5j which compares the magnetic 
structure on the final day of periods A (top two panels) and C (bottom two 
panels), representing minimum and maximum respectively. The upper panel in 
each pair shows flux rope points overlayed on the photospheric magnetic field 
(with differently coloured clusters corresponding to different flux ropes), while 
the lower panels show the distribution of current helicity (density) a = jo • Bo / 
over the spherical surface at height 12 Mm, in the low corona. Flux ropes arc 
seen to form at locations where current helicity is concentrated above photo- 
spheric PILs, which is a key feature of this model. It is therefore natural that at 
maximum — when there are more PILs — there should be more flux ropes. Indeed, 
both the total PIL length and the integral of |a| over this spherical surface show 
increases of about a factor 1.5 from minimum to maximum, comparable to the 
number of flux ropes. 

We remark that the precise origin of the magnetic helicity that becomes 
concentrated above PILs is stil l an open question, and ne ed not be the same 
for PILs everywhere on the Sun. lYeates and Mackavl (l2009al ) found that sheared 
magnetic fields overlying PILs in this model arise purely from differential rotation 
when the PILs are at high latitudes, but over low-latitude PILs shear may also 
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Figure 5. Comparison of low-coronal magnetic field on final day of periods A and C. The top 
panels for each period show the photospheric radial magnetic field (solid contours/ white for 
positive, dashed contours/gray for negative) and locations of fiux ropes on this day (coloured 
points). The bottom, red shaded, panels show the distributions of current helicity a at height 
12 Mm (white indicates positive, black negative, saturation level ±20 X 10~*m~-'^), with 
photospheric PILs in blue. 
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arise from the emergence of twisted magnetic fields in active regions, and from 
the interaction between regions of strong magnetic field. 

A consequence of the dependence of flux rope locations on PILs is that the 
latitude distribution of flux rope points varies over the solar cycle. This is shown 
by the thick histograms in Figure [HI At minimum (periods A and F) there 
is a bimodal distribution at active latitudes, whereas in periods B to E the 
distribution of flux ropes extends both across the equator and to higher latitudes. 
High latitude (above 60°) flux ropes are particularly prevalent in period B. This 
is due to the poleward migration of the polar crowns, prior to the polar field 
reversal. 



3.3. Flux Rope Ejections 

To detect flux rope eje ctions in each simulation ru n, we have again used the 
procedure described by lYeates and Mackav selecting those flux rope 



points with -uor > 0.5kms~^ in the magneto-frictional code, and clustering in 
both space and time into separate ejection events. With the additional factor _Bq 
in the detection of flux rope points (Section 13. 2|) , we found that performance of 
the ejection detecting algorithm was optimized if the parameters were slightly 
adjusted: we now require that clusters are separated by at least 5 days in time, 
and have at least 8 points. 

To look at how the simulated ejection rate varies over the solar cycle, we 
count the number of ejections in the last 100 days of each simulation period. 
This varies from a minimum of 11 in period A up to a maximum of 94 in period 
B. Not e that this ejection rate for p eriod B is lower than those in the simulation 



runs of lYeates and Mackavl (l2009bf ) which had uniform signs of bipole twist of 
magnitude |/3| = 0.2 or 0.4 in each hemisphere. In the present simulation, nearby 
bipolcs with opposite signs of P tend to mutually "cancel" some of their helicity. 
Figures Ulja) and (b) show how this simulated ejection rate varies over the solar 
cycle, respectively normalised (by the value for period A) and un-normalised. In 
this case the error bars show estimated errors i n the number of ejec t ions re sulting 



from the automated detection procedure see lYeates and Mackavll2009b[ ). 

This eightfold increase in the rate of ejections is similar to that of the total 
magnetic energy in the corona, or that of the bipole emergence rate. We now 
make two comparisons. Firstly, the cycle increase in ejection rate is larger than 
that in the number of flux rope points present at any one time (Section 13. 2p . 
Therefore a higher proportion of flux ropes lose equilibrium in a given interval 
at maximum than at minimum, so fltxx ropes at maximum have shorter lifetimes 
on average. Secondly, the increase in ejection rate is only about half that in free 
magnetic energy (Section [XT]) . Hence the total free energy, as plotted in Figure 
[IJb), does not determine directly the number of flux rope ejections occurring, 
but rather the local structure and evolution of the magnetic field must be taken 
into account. In particular, the similar increases in ejection and bipole emergence 
rates suggest that newly-emerging active regions are important. 

We can see that the increased rate of ejections is somehow related to newly- 
emerged regions by inspecting the latitude distribution of ejection points (the 
filled histograms in Figure E]). In Periods B to E, the distribution of ejection 
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Figure 6. Histograms showing latitude distribution of flux rope points (thick lines) and 
ejection points (filled), for final 100 days of each simulation period. Bin sizes arc 4° and 
5° respectively. 



points is more bimodal than that of flux rope points, and has fewer points 
at higher latitudes, with only a few ejections on the polar crowns. Indeed, a 
similar bimodal latitude distribut ion is found in observation s of disappearing fil- 
amen ts ( Poioga and Huaii3 . 2003 ). of prominence eruption s (iGopalswamy et~al. 



20031). and of EUV sour c e regi ons associated with CMEs (jPlunkett et all 12001 



Cremades and Bothmeii 12004 ). 



There are several mechanisms by which increased flux emergence causes lo- 
cally increased propensity for flux rope ejection. Firstly, increased emergence 
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injects a greater amount of magnetic energy and helicity, leading to more rapid 
growth of flux ropes and hence earlier loss of equilibrium. This is consistent with 
the finding of similar relative increases in both ejection rate and total magnetic 
energy. Secondly, active regions have relatively more magnetic energy than the 
overlying field, so flux ropes in or around active regions are more likely to lose 
equilibrium. Finally, there is the possiblity of destabilisation of existing flux 
ropes by the nearby emergence of new active regions. 



4. Observed CME Rate 



How does the rate of flux rope ejections produced by the model compare to the 
observed CME rate? First, note that care should be take n with observed CME 
rates , because these depend both on instrument sensitivity ( Cremades and St. Cvij .[ 



2007p and on t he pr ocessing of the data, including the definition of a CME {e.g., 



Yashiro et aLl . 120081 ). 



These limitations notwithstanding, we have determined two observed rates 
over cycle 23, based on the published CDAW and CACTus catalo gues, both using 
data from LASCO (Large Angle and Spectrometric Coronagraph; Brueckner et al\.t 
19951) . The CDAW (Coordinated Data Analysis Workshops) catalogue ( Yashiro et aZ.I .[ 



2004L available at http://cdaw.gsfc.nasa.gov/CMEIist/) is the standard manually 



compiled list of LASCO CMEs; we filter out CMEs with apparent width less than 
15° or greater than 270°. The resulting number of CMEs per day is shown by the 
red curves in Figures SJ^a) and (b). The blue curve in Figure Hl^b) shows the rate 
from the alternative CACTus catalogue (http://sidc.be/cactus), which is based 
on the same LASCO observation s but uses an automated technique to detect ra- 
dial motion in height-time maps ( Robbrecht and Berghmai5 2004). We omit so- 



called "marginal cases" in the CACTus data (see http://sidc.oma. be/cactus/scan/).| 

The error bars on these observed rates in Figure H] are d etermined tak i ng int o 



account instrument data gaps, following the method of ISt. Cvr et al. (I2OOOI) 



Namely, the lower bar is the actual observed number of CMEs, the data point 
itself has additional CMEs at this same rate to fill any data gaps, and the upper 
bar fills the data gaps with CMEs at the maximum rate observed over any 
24-hour period. 

Figure SKa) shows an appealing agreement between the relative increases of 
simulated and observed (CDAW) ejection rates over the cycle. However, Figure 
SKb) shows that the actual simulated rate is much lower than the observed 
rates from either catalogue (which agree quite well). The ratio of simulated to 
observed ejection rates remains roughly constant, at about 0.2, over the cycle. 
The implications of th is shortfall fo r possi ble mechanisms of CME initiation are 



considered elsewhere (Yeates et al. 20101). Note that this ratio of 0.2 is lower 



than that found by lYeates and Mackavl ( 2009lJ ) in their optimum case; partly 



because they omitted "poor events" in the CDAW catalogue, and partly because 
the bipoles in this paper have been given random twists, rather than all having 
the same value of /? a s in the earlier simulations. T his results in lower net helicity, 
which was shown bv lYeates and Mackavl (|2009bl ) to result in fewer eruptions. 
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5. Conclusions 

As a step toward removing tiie potential-field assumption in modelling the global 
coronal magnetic field, we have investigated a quasi-static model allowing for 
electric currents. This model is driven by photosphcric magnetic observations 
and follows the build-up of magnetic helicity as magnetic flux emerges and is 
transported by surface motions. In this paper, we simulate six distinct periods 
during solar cycle 23, so as to study how the magnetic field structure predicted 
by this model varies over the solar cycle. 

A key feature of this model is the transport of magnetic helicity and its 
concentration into twisted magnetic flux ropes above photosphcric PILs. The 
solar cycle variation of these flux ropes in the model may be characterised as 
follows: 

1. The number of flux ropes present at any one time doubles between minimum 
and maximum of the solar cycle. This follows the total length of photosphcric 
PILs. 

2. The rate of flux rope ejections — or losses of equilibrium — increases by a factor 
of eight between minimum and maximum. This is lower than the relative 
increase in free magnetic energy, but similar to the relative increases in total 
magnetic energy and total parallel current, both of which follow the number 
of emerging active regions. 

The difference between the relative increases in ejection rate and in the number 
of flux ropes implies that, at maximum, the flux ropes have shorter lifetimes 
before losing equilibrium. This is a consequence of the higher rate of active 
region emergence. 

Although there are fewer flux rope ejections at minimum, it is not the case 
that the corona in our model is everywhere close to potential. Even in the 
absence of many emerging bipoles, shearing of pre-existing coronal field by 
differential rotation generates current on a large scale at all latitudes. A weak 
background field originating in earlier decayed active regions is present even 
in the recent extended minimum (modelled by our period F in 2008). Our 
model assumes that differential rotation is equally effective at shearing this 
weaker field, so that, for example, we find systematic sheared arcades on the 
high-latitude polar crowns. Interestingly, although sheared arcades are formed, 
the low reconnecting flux over these PILs in 2008 prevents detached helical 
flux ropes from forming at many locations, unlike at solar maximum. The true 
nature of the magnetic field above the polar crown PILs remains uncertain. Our 
simulations predict that differential rotation will develop positive and negative 
helicity over such east-west PILs in the northern and southern hemispheres 
respectively. However, the (few) existing observations of the magnetic field in 
polar crown prominences imply the opposi te sign of helicity in each hemisphere 



(jLerov. Bommier and Sahal-BrechotL 119831 ). In our quasi-static model, the cor 



rect sign of helicity is onl y obtained at such PIL locations if additional emergence 
of axial flux is included ( van Ballegooiien. Priest, and Mackavl [200(1 . This was 



not included in the present paper because it lacks observational justification. A 
concerted effort to obtain further observations of the magnetic field structure 
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at high latitudes is needed in order to resolve this issue and hence test whether 
additional sources of helicity are required in addition to those in our existing 
model. 

An important question that may be addressed by a time-dependent non- 
potential model of the coronal magnetic field is the initiation of CMEs. Our 
model allows us to test whether the loss of equilibrium of flux ropes formed 
by quasi-static shearing and flux cancellation is sufficient to account for the 
observed CME rate. For the present form of the model the answer is negative: 
the ratio of simulated to observed ejection rate (from CDAW) remains about 0.2 
over the whole solar cycle. Even in our earlier (less realistic) simulations where 
all emerging bipoles were given the same sign of helicity in each hemisphere, the 
ratio increases only to about 0.25. This suggests strongly that the formation of 
flux ropes by quasi-static shearing of the large-scale magnetic field is incapable 
of initiating all CMEs. We should point out, however, that the CDAW rate for 
period F (in 2008) may be too high (as discussed in Section . in which case 
the model may account for a greater proportion of CMEs in the recent very 
quiet period. Nevertheless, it is clear that the present model lacks the detail to 
adequately describe initiation of all CMEs. 

More detai led comparison wi th observed low-coronal CME source regions 
(described by lYeates et al. , 2010[) reveals that many of the additional observed 



CMEs are located in active regions, with individual active regions frequently 
producing multiple CMEs in the same day. Our global quasi-static simula- 
tions cannot produce such dynamic events with the present form of input data 
(updated once per Carrington rotation). The quasi-static approach could be 
developed in future to incorporate driving from photospheric magnetograms on 
much shorter spatial and time scales, in order to test whether the energetic 
restructuring of the corona following flux emergence may be adequately described 
by a quasi-static model, or whether fully dynamic MHD modelling is required. 
From the global perspective, this more detailed modelling of active regions should 
allow their helicity content, and hence the global transport of helicity, to be 
better constrained by observations. This has important implications both for 
models of the coronal magnetic structure and, ultimately, for the origin of the 



Sun's magnetic field in the solar dynamo (jSeehafer et all l2003r ) 
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